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ANALXTIC SOLUTIONS TO THE IGNITION KINETICS 


OF THE HYDROGEN-OXYGEN REACTION* 
by Richard S. Brokaw 
Lewis Research Center 


SUMMARY 

Solutions have been obtained by assuming isothermal conditions and negli- 
gible depletion of reactants during the induction period for the following 
scheme of reactions : 

k i 

OH + H 2 — ►HgO + H 
k 2 

H + 0 2 — ^OH + 0 
k 3 

0 + H 2 *-0H + H 

H + 0 2 + M— ^-H0 2 + M 
k ll 

H0 2 + H 2 H 2 0 2 + H 

Three types of solutions have been deduced corresponding to kgO^ > 2k 2 (a 
low-temperature - high-pressure region of long ignition delays), 2 k 2 > k 0 C^ (a 
high -temperature - low-pressure region of short delays ), and 2 k 2 = ^Q c y[ ("the 
boundary between the two regions). Experimental ignition delays in both long- 
and short-delay regions are adequately explained, and ignition delays computed 
by numerical integration of the rate equations are reproduced to within a few 
percent. Finally, ignition delays for the Von Neumann spike condition in 
Chapman -Jouguet detonations in the neighborhood of the lean limit of deton- 
ability in air have been calculated* These delays decrease by two orders of 
magnitude as the hydrogen mole fraction increases from 0.135 to 0.145. This is 
very close to the experimental limit (hydrogen mole fraction, 0.14 to 0.15) and 
thus supports the notion that the condition 2 k 2 > kgc^ must be satisfied for 
a stable detonation. 


INTRODUCTION 

When hydrogen and air are instantaneously mixed at high temperature (as in 


*This report was 
posium on Combustion, 


also presented as a paper at the Tenth International Sym- 
Cambridge, London, August 17 to 21, 1964. 



idealized analyses pertinent to hypersonic ramjets, refs. 1 and 2), or when a 
mixture containing hydrogen and oxygen is subjected to a pulse of high tempera- 
ture and pressure (as in shock tube ignition studies, refs. 3 and 4), the fol- 
lowing events ensue: First, there is a very short period during which small 

concentrations of atoms and radicals - H, OH, and 0 - are produced. The mech- 
anism of this initiation is uncertain and may vary under different circum- 
stances. There follows the hulk of the induction period, during which the 
chain reaction is propagated and branched through the elementary reactions 

k l 

OH + H 2 *-H 2 0 + H I 

H + Q 2 . k g».OH + 0 II 


k 3 

0 + H 2 — *-0H + H III 

The first reaction is not strongly exothermic (i.e., heat of reaction, 

15 kcal/mole at 1000° K) , while the second and third reactions are endothermic 
by 16 and 2 kilocalories, respectively. Thus, during most of the induction 
period, the atom and radical concentrations increase exponentially under es- 
sentially isothermal conditions. Very late in the induction period atom and 
radical concentrations become so large that self -heating occurs, which is 
largely due to the highly exothermic three-body recombination processes. The 
temperature rises precipitously, and the induction period is at an end. 

At the lower temperatures and higher pressures the further reactions 


k fi 

H + 0 2 + M — ^-H0 2 + M 
k ll 

H0 2 + H 2 *-H 2 0 2 + H 


VI 

XI 


must be considered as veil. These five reactions, together with terms describ- 
ing the destruction of chain carriers after diffusion to the vails, suffice to 
explain qualitatively the ignition behavior of the hydrogen-oxygen system, in- 
cluding the existence of the first, second, and third explosion limits. 


The differential equations describing the chemical kinetics of reactions 
I, II, and III have been integrated numerically for hypersonic ramjet (refs. 

1 and 2) and shock -tube (ref. 4) conditions. This report presents explicit 
analytic solutions that are excellent approximations to the numerical re- 
sults; furthermore, reactions VI and XI have been included. Similar analytic 
solutions for reactions I to III and VI have been obtained by Kondratiev 
(ref. 5); Asaba, Gardiner, and Stubbeman (ref. 6) used the solution for re- 
actions I to III to analyze results of their shock -tube experiments. The ex- 
plicit analytic solutions developed in this report are valuable in that they 
permit one to determine easily the effects of variations in specific rate 
constants, temperature, pressure, initiation mechanism, etc., on ignition de- 
lays and the like* Some illustrative applications are presented and dis- 
cussed* 
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THEORETICAL CONSIDERATIONS 

The differential equations governing the growth of radical concentrations 
during the induction period are as follows: 

dc 0H _ 

dt k l c H2 C 0H + k 2 c 02 C H + k 3 c H2°0 + 1 0H (!) 

dcjj 

dt~ = k l c H 2 c OH " k 2 c 0 2 C H + k 3 c H 2 c 0 " k 6 c 0 2 c M c H + k ll c H 2 c H0 2 + % ( 2 ) 

d c 0 

dt“ = k 2 c 0 2 c H “ k 3 c H 2 c 0 + 1 0 ( 3 ) 

dc H0 2 

“dt = k 6 c 0 2 c M c H " k ll c H 2 c H0 2 + i H0 2 ( 4 ) 

Here, the k's are the specific reaction rate constants, the c’s are con- 
centrations^ and the i's are rates at which the radicals are generated spon- 
taneously, either in the gas phase or on the surfaces* (it is also possible to 
include destruction of radicals at the walls hy including terms like T v CoHJ 
etc., in eq_s. (l) to (4) with the coefficient , a function of diffusivity, 
pressure, geometry, and surface activity. ) 

During most of the induction period the kinetics of the elementary reac- 
tions I to III, VI, and XI can be represented as pseudo-first-order reactions; 
that is, they are all first order in the atom and radical concentrations - the 
only concentrations that are changing appreciably. 

The initiation rates i()Hj etc., which are zero order in the radical con- 
centrations, can be eliminated from equations (l) to (4) by means of the sub- 
stitutions (ref. 7) 



These initiation rates are presumably very small compared with the rates of re- 
actions I to III, so that the correction terms of equations (5) to (8) are 


3 



negligible except in the very initial portion of the induction period. Con- 
sequently in what follows these corrections will he neglected. They should he 
considered, however, in any detailed analysis of the kinetics of initiation. 
Inclusion of these terms poses no mathematical problem; they simply increase 
algebraic complexity. 

The system of differential equations (eqs. (l) to (4)) may be solved by 
the method outlined by Frost and Pearson (ref. 8) if it is assumed that the 
concentrations of reactants - Hg, O 2 and M - are constant and neglecting i-0H> 
ijj, i Q, and igo 2 (or applying eqs. (5) to (8)). A set of particular solutions 
is 

c ± = A i e At i = OH, H, 0, H0 2 (9) 

where is a constant with dimensions of concentration and A is a constant 

with dimensions of reciprocal time. Substituting equations (9) into equations 
(l) to (4) yields the relations 


" ( k l c H 2 + 

+ k 2 c 0 2 A H 

+ k 3 c H 2 A 0 

= 0 

(10) 

k l c H 2 A 0H " 

( k 2 c 0 2 + k 6 c 0 2 c M + A) A H 

+ k 3 c H 2 A 0 

+ k ll c H 2 A H0 2 = 0 

(11) 


k 2 c 0 2 A H 

" ( k 3 c H 2 + A)Aq 

= 0 

(12) 


k 6 c 0 2 c M A H 

- 

( k ll c H 2 + A ^ A H0 2 = 0 

(13) 

A nontrivial solution requires that the 

determinant of the coefficients 

of the 


Aj_ be zero (ref. 8). This leads to the quartic equation 
+ |( k l + k 3 + k ll) c H 2 + ( k 2 + k 6 c M) c 0 2 J^ 3 

+ |k 1 k 3 c| 2 + (kq + k 3) k 6 c H 2 c 0 2 c M + ( k l c H 2 + k 2 c 0 2 + k 3 c H 2 ) k ll c H 2 J ^ 

+ k i k 3 c| 2 [(k 6 c M - 2k 2 )c 02 + kqqc-gjj A - 2kqk 2 k 3 kqqcg 2 c 02 =0 (14) 

From Decartes' rule of signs, it is found that this equation has one positive 
root, one negative root, and two roots that are either negative or conjugate 
complex. The positive root is desired since it governs the exponential growth 
of radical concentrations during the induction period. The other roots are 
important only in the initiation phases of the reaction - their effect on the 
radical concentrations dies out (more or less quickly depending on their 
absolute magnitude). The rate constant kqq is very small in comparison with 
kq, k 2 , k 3 , and k g c M . Therefore the additive terms involving kqq in the 
coefficients of A, A 2 , and A^ may he neglected. 

Approximations to the positive root of equation (14) will he derived for 
three conditions: kgc^ > 2k 2 , k g c M = 2kg, and k g c M < 2kg. 
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kg c M > 2kg. - This regime corresponds to the conditions above the third 

explosion limit. This is a region of long ignition lags, with delays of sev- 
eral seconds at 560° to 600° C and atmospheric pressure. Correspondingly, the 
positive A is very small, and it is found that A 2 , A 3 , and A^ may he set 
equal to zero in equation (14). The result is 


A = 


2k 2 k ll c H2 
k 6 c M " 2k 2 


(15) 


This expression is a very good approximation indeed, except very close to the 
second limit condition, where the denominator approaches zero. Note that at 
high pressure A — ^ 2 ^2^11 C H /^6 C M* is independent of pressure. 


kg c M = 2k2 (Second limit condition). - Above the junction of the second 

and third limits this condition marks the boundary between regions of long and 
short ignition delay. Here the coefficient of A in equation (14) is zero 
(neglecting the term in k^). Again, A 3 and A 4 maybe neglected; the re- 
sult is 

2k 1 k 2 k 3 k 11 c| c 0c I 1 / 2 


A = 


h 


2 >J 2 

^l k 3 c H2 + ( k l + k 3) k 6 c 02 c M 


- J 

l 


(16) 


2kg > kgcjq. - This corresponds to the important region of short ignition 

delays at high temperatures and modest pressures. Under these conditions the 
coefficient of A in equation (14) reduces to the cubic equation 


k l k 3 1 

( 

[f k 2 - 

k 6 c M > 

|c h 2 c o 2 

k,k c 2 + (k n + k )k c c ■ 
1 3 H 2 1 3' 6 H 2 C>2 

"M + 

( k -j_ + 

k 3 )c H 2 + <k 2 + k 6°M> C 0 2 

A + A 2 


which was first derived by Kondratiev (ref. 5). Note that, if A is small, 
the terms involving A and A 2 in the denominator are unimportant. If, in 
addition, kq, k2, and k3 are much larger than kgCjyj, A approaches ^^cq^. 
This provides the basis of the ignition delay correlation used by Scott and 
Kinsey (ref. 3). 

Equation (17) has a further interesting aspect. Suppose one has an ap- 
proximate value of A, designated Aq. The right side of equation (17) can 
thus be used to generate a second approximation, A2* A positive root is 
needed. Hence, if Aq is too small, Ag is too large; if Aq is too large, 
A2 is too small. Thus, the correct solution is bracketed, and the approxima- 
tion A = (AqAg) 1 / 2 should be in error by less than half the error of the 
worst of the two estimates. i 


iNote that, if the error in A2 is less than twice the error in Aq, the 
geometric mean value (call it A3) is clearly superior to Aq. Thus, A3 can 
be substituted into eq. (17) to obtain an improved approximation, Aq; and 
Ag = (A3A4) 1 / 2 wiii be an improvement over A3 = (AqAg) 1 / 2 . This provides 
an iterative procedure that wiil converge on the true vaiue of A. 
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An approximate solution can "be obtained by invoking the steady-state ap- 
proximation in OH (by setting eq. (l) equal to zero). Thus, equation (l) can 
be used to eliminate the OH concentration from equations (2) and (3), which can 
then be solved with the result 


An 


k 3 cr H2 + k 6 co 2 c M I 


4k 3 (2k2 - k6CM)cH 2 co 2 1 1 / 2 


1 + -7 - lb 

( k 3 c H 2 + k 6 c 0g C Mj 


) " 


A 


= (A^g) 1 / 2 


k- L k 3 ( 2k 2 - k 6 c M ) c | 2 c 0g 


{P 


4k 3( 2k 2 ~ k 6 c M) c Hp c O P "| 1//2 


1 ( k 3 c H 2 +k 6 c 0 2 c M) ? ' 


(18) 


This result, in turn, may be substituted into equation (17) to obtain A 2 , so 
that the final approximation to A is 


k 3 c H 2 _ k 6 c 0 2 c M 

V 2k 2 c 0 2 k 3 c H2 + k 6 CQ 2 c M + ^ k l c H 2 + k 2 c 0' 



2 co 2r 1 

)2 J l| 

■ k^CM) ctt cn_ Hi/ 2 


4k3( 2k 2 " k 6 c M) c Hg C Og 


k 3 c H, C. 


e^Og^M 


n 


4 1 


(19) 


It is possible, of course, to obtain other first approximations A^_ by- 
assuming steady states in H or 0 rather than in OH. The following table shows 
some numerical comparisons with these other techniques calculated for a mixture 


Steady state 
in - 

*'_i 

Temperature , °K 

sec 

1100 

1900 



Pressure, atm 



0.2307 

0.4873 

OH 


2.67X10 4 

3.70X10 5 


A 2 

2.39 

2.02 


<V 2 > 1/2 

2.526 

2.734 

H 


10.1X10 4 

4.58X10 5 


*2 

.685 

1.616 


(A-jAg) 1 / 2 

2.630 

2.720 

0 


8.76X10 4 

16.8X10 5 


*2 

.793 

.360 


(A!A 2 ) 1/2 

2.638 

2.460 


True value 

2.503X10 4 

2.727X10 5 


of 5 percent hydrogen in air; the rate 
constants are from reference 4, and 
kg has been set equal to zero. The 
assumption of a steady state in OH 
gives clearly superior values of A^, 
although the final approximations 

are remarkably good in all 
cases. The OH concentration is smal- 
lest, so it seems reasonable that the 
steady-state approximation should ap- 
ply best to this radical. (The true 
values of A were obtained by numer- 
ical solution of eq. (17).) 

Once the value of A has been ob- 
tained, relations among AqHj Ah, Aq, 
and Aho 2 follow from equations (10) 
to (13). 

It remains merely to discuss the 
boundary conditions defining the be- 
ginning and the end of the ignition 
delay period. For simplicity it will 
be assumed that the initiation phase 
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of the induction period is very short and may he neglected in comparison with 
the total span of the delay before ignition. (As indicated heretofore, it is 
possible to obtain analytic solutions without this assumption, but at the ex- 
pense of algebraic complexity. Further a detailed knowledge of the initiation 
kinetics would be required*) 

Assume that the initiation process produces radicals by some bimolecular 
process : 

H = (dT7 t =o K ( 2 °) 

where c^ represents the total concentration of stable species. Examples of 
reactions of this sort could include 

H 2 + 0 2 ^2 OH 

H 2 + 0 2 ►H + H0 2 

H 2 + M >-2H + M 

and others as well* It is reasonable to assume that, initially, the rate of 
formation of radicals is equal to the rate of consumption In reactions such as 
I, II, III, etc*, (which are first order in the radical, but second order over- 
all) : 

(dt ) t=0 “ c M A i ( 21 ) 

where Ai is an effective initial radical concentration (see eq. (9)). Thus, 
in general, from equations (20) and (21), 

A^ = const cj 4 = aj_p (22) 

where p is pressure and aj_ is an appropriate constant of proportionality. 
Hence, the radical concentration as a function of time is 

c-j_ = aj_pe^ (23) 

as long as initiation is in the gas phase and requires only a small fraction of 
ignition delay* 

The end of the Induction period can be defined in a number of ways that 
are qualitatively equivalent* For example, Schott and Kinsey (ref* 3) define 
the end of the delay period as the time at which the concentration of some 
species reaches a characteristic level (in their case, cqjj = 10" ^ mole/liter). 
Thus, from equation (23), 

At = B - In p (24) 

where t Is the ignition delay and 33 is a constant* The same result is ob- 
tained if one characterizes the end of the delay period by a critical rate of 
temperature rise due to three-body radical recombination 

C M C p (if) crlt * c i c j c M K c MP 2e2AT (25) 
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Ill 


which again leads to equation (24). Again, if one assumes a critical rate of 
temperature rise due to a bimolecular process such as reaction I, equation (24) 
follows once more. Thus, a variety of reasonable assumptions about the initia- 
tion and the termination of the induction period all lead to expressions of the 
form of equation (24). 


ILLUSTRATIVE APPLICATIONS 

In this section, the foregoing theoretical results are applied to cal- 
culate ignition delays for hydrogen-air systems under a variety of conditions 
of practical interest, including (l) the low-temperature region of moderate 
pressures, where ignition delays of one to several seconds are observed, 

(2) the higher-temperature low-to-moderate -pressure regime, where delays range 
from several microseconds to several milliseconds, and finally (3) the high- 
temperature - high-pressure boundary between these regimes, which seems to be 
of importance in determining limits of detonability. 


Region of Long Delays 

Ignition delays for hydrogen in air were measured some three decades ago 
by Dixon (ref. 9), who used an apparatus in which a heated stream of hydrogen 
was injected into a larger concentric stream of heated air. Interpolated de- 
lays at 600° C are shown in figure 1. 

In this regime, the characteristic parameter A (used in eq. (24) to com- 
pute the ignition delay t) is given by equation (15). Prom figure 1 it is 
clear that, at a pressure close to l/2 atmosphere, the ignition delays become 
very short, as kgcjj approaches 2k 2 . This provides a key for determining the 
effective hydrogen concentration (which is otherwise indeterminate in Dixon’s 
experiment), because kg is a function of the nature of the third body, M. 

It is assumed here that 

kg = 3.27xl0 15 (x H2 + 0.35 x Q g + 0.43 x N +0.2 x Ar 

+ 14.3 + • • . )T -1 ' 92 liters 2 /(mole) 2 (sec) (26) 

where the x’s are the mole fractions of various third bodies. Equation (26) 
is consistent with Clyne’s room temperature measurement, with argon as a third 
body (ref. 10), and Baldwin's estimates at 540° C (ref. 10). The relative 
third-body efficiencies are taken from reference 11. Using kg from equa- 
tion (26) and kg from reference 4 gives an effective mole fraction of 
hydrogen of 0.22. This procedure may seem somewhat arbitrary, but it is 
comforting that the hydrogen mole fraction did not turn out to be negative 
(limit pressure > 0.66 atm) or greater than unity (limit pressure < 0.27 atm) i 

The curve in figure 1, computed by assuming k^ = 5.2xl0 2 liters per mole 
per second and B = 3.8, provides a very good fit of the experimental data. 

This value of k^ is greater than that suggested by Voevodsky (ref. 12) by a 
factor of 4.5. The cause of this discrepancy is uncertain: Voevodsky perhaps 

deduced his kqq by assuming somewhat different values of kg and kg; but on 
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Ignition delay, sec 


the other hand our kinetic scheme is probably oversimplified. Furthermore, the 
effective hydrogen concentration may vary with pressure. 


Other observations by Dixon (ref. 9) are also easily understood from these 

theoretical relations. Thus, he reported 
ignition delays in oxygen slightly short- 
er than those in air. This is not sur- 
prising since the oxygen concentration 
does not appear explicitly in equa- 
tion (15); its only effect is to cause a 
slight decrease in kg (see eq. (26)). 

On the other hand, Dixon (ref. 9) re- 
ported markedly longer delays in the 
presence of water vapor. This is due to 
the fact that water is a very efficient 
third body , so that kg is considerably 
increased. 


The effect of hydrogen concentration 
on the ignition delay has been determined 
by Anagnostou, Brokaw, and Butler 
(ref. 13); some of their results near 
600° C are shown in figure 2. (They 
found some variation according to the 
mode of operation of their apparatus; the 
shortest delays at this temperature were 
roughly half as large as those shown in 
has been calculated from equations (15) and (24) 
fit Dixon *s data in figure 1. Thus, the ex- 
perimental composition dependence is 
rather well predicted: t cc p“0* 



fig. 2.) The curve in figure 2 
with the same constants used to 



l 

.05 


i 

.10 .15 .20 .25 .30 

Mole fraction of hydrogen 


.40 .50 


Figure 2. - Effect of hydrogen concentration on hydrogen-air 
ignition delays at 594° C and 1 atmosphere (ref. 13). 


It seems then that the analytic 
solutions presented in this report 
describe the behavior of the 
hydrogen-oxygen reaction in the re- 
gion of long delays rather well. 

Region of Short Delays 

Ignition delays for hydrogen-air 
have been calculated numerically by 
Momtchiloff, Taback, and Buswell 
(ref. l) and Belles and Lauver 
(ref. 4) with different assumptions 
as to the mechanism of initiation. 
Momtchiloff, Taback, and Buswell use 
Schott and Kinsey's criterion for 
ignition (c qh = 10“ 6 mole/liter. 


ref. 3); Belles and Lauver (ref. 4) give constants for equation (9) so that 
Schott and Kinsey's criterion is easily applied to their calculations as well. 
(Both papers neglect reaction VI. ) 
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The analysis of reference 1, presumably pertinent to hypersonic ramjets, 
assumes equilibrium concentrations of hydrogen and oxygen atoms in the hydrogen 
and air streams prior to instantaneous mixing. The hydrogen atoms are expected 
to be more important because of their higher equilibrium concentration. Thus, 
the growth of OH radical concentration with time is 


. _ a oh 

C 0H = c H,eq 5 ^ 


e At 


(27) 


Using equations (10) and (11)., neglecting reactions VI and XI, gives 


A 0H 2k 2 cp 2 + A 
a E ~ 2k l c H2 + A 


Thus, the induction times computed numerically in reference 1 can he approxi- 
mated from equations (19), (27), and (28). numerical and approximate analytic 
results are compared in the following table. 


Temperature, 

°K 

Pressure, 

atm 

Mole fraction 
of hydrogen 

Ignition delay, psec 

Approximate 

(a) 

Numerical 

973 

1.0 

0.296 

240 

b 260 

1667 

1.0 

.296 

5.3 

b 5.1 

1100 

.2307 

.05 

866 

c 866 

1900 

.4873 

.05 

33.6 

c 32.6 


a This report. 
b Ref. 1. 

C Ref. 4. 


Belles and Lauver 
(ref. 4), on the other 
hand, assume initiation 
by the reaction 


H 2 + 0 2 



20H 


If it is assumed that, 
initially, the rate of 
production of OH radicals 
in this reaction equals 
the rate of their con- 
sumption in reaction I, 


aoh 


2k 0H c 0 2 


(29) 


so that the ignition delay is 


At = In 



k-]XlO _6 \ 
2k 0H c 0 2 / 


(30) 


Numerical (ref. 4) and analytic results from equation (30) are compared in 
the previous table. It is clear that the analytic solutions are very good ap- 
proximations to the numerical results, both for very lean and stoichiometric 
conditions, and for two rather different assumptions as to the initial source 
of radicals. 

By studying stationary hydrogen-air flames in a flow system, Foure 
(ref. 14) has obtained information on ignition delays in the region of short 
delays. Foure heated flowing homogeneous fuel-air mixtures to a temperature 
somewhat below the point of spontaneous ignition. The mixture, traveling at a 
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high subsonic velocity, then passed through a diffuser into a reaction tube. 

The velocity decrease to about 40 meters per second was accompanied by in- 
creases in temperature and pressure that rendered the mixtures spontaneously 
inflammable. Temperature and composition profiles were obtained by probing the 
extended reaction zone formed some distance downstream of the diffuser* The 
distance from the diffuser exit to the point of appreciable temperature rise 
provides an indication of ignition delay times; however, the choice of the dif- 
fuser exit as a zero point is somewhat arbitrary since reaction may actually 
start within the diffuser. Foure f s results are shown in figure 3. A milli- 


second scale has been added based on the 
reaction tube was 40 meters per second. 



Figure 3. - Effect of hydrogen concentration on hydrogen- 
air ignition delays for 658±2°C and 1 atmosphere (ref. 14). 


cur at 912° or 913° K. Thus, the 
perature limit within 3° K. 


assumption that the velocity in the 

Ignition delays computed from 
equations (19) and (24) are shown by 
the curve in figure 3. The constant 
B has been taken as 25 so as to fit 
the experimental data at the leaner 
compositions. (This agrees very well 
with eq. (30), which predicts 
B = 25. 6 at 913° K. ) The deviation of 
the two richest mole fractions sug- 
gests that there may have been consid- 
erable reaction in the diffuser under 
these conditions. 

Foure was unable to observe any 
temperature rise in his ignition tube 
below 910° K. Since the tube was 1 
meter long and the flow velocity was 
about 40 meters per second, ignition 
delay must have been in excess of 25 
milliseconds. If B = 25 is again 
assumed, an ignition delay of 54 mil- 
liseconds is computed. At this condi- 
tion, the delay is extremely sensitive 
to the temperature. If the tempera- 
ture is raised but 1 degree the delay 
decreases to 40 milliseconds , and 
hence a 25-millisecond delay would oc- 
Lt analysis predicts Foure T s lower tem- 


A further point merits some discussion. The ignition delays in the re- 
gion of long delays were fit assuming B ^ 3. 8 at 600° K, whereas the short 
delay data require B ^ 25 at a temperature some 60° higher. From equations 
(20) to (25), it appears that 


B = - In 



+ 2 In p + const 


so that a small value of B indicates a large initiation rate and vice versa. 
Thus, it appears that the initiation mechanisms in the long- and short-delay 
regions are distinctly different, and the long delay initiation rate is larger 
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by perhaps nine orders of magnitude (although this changes the delays by less 
than one order of magnitude). There is evidence (ref. 15) that hydrogen perox- 
ide plays an important role in the initiation of the hydrogen-oxygen reaction 
in the neighborhood of the second limit. Hydrogen peroxide is produced in re- 
action XI, which is important principally in the long-delay regime. 

Finally, it should be noted that this approximate solution for the region 
of short ignition delays adequately represents shock-tube ignition data on the 
hydrogen- oxygen system over a very wide range of mixture ratio (ref. 16). 


Boundary Between Long- and Short-Delay Regions 
(and Limits of Detonability) 

A stable detonation wave involves the coupling of a shock front with a 
rapid chemical heat release, which, in turn, sustains the shock. Thus, the 
heat release must occur rather close to the shock in a stable detonation. 

Belles (ref. 17) has proposed that, in the case of the hydrogen- oxygen re- 
action, the limits of detonability may be predicted by requiring that in a sta- 
ble detonation the second limit condition 2k2 > kgcjyi be satisfied in the 

Von Neumann spike accompanying the detonation front. (These conditions of 
temperature and pressure are calculable from the Chapman- Jouguet detonation 
Mach number, and the adiabatic shock relations for the nonreacting gas at the 
same Mach number. ) Using the second limit criterion. Belles predicted lean and 
rich limits of detonability for hydrogen-oxygen and hydrogen-air in good agree- 
ment with experiment. 

Belles applied the second limit criterion under conditions of temperature 
(~1300° K) and pressure (~20 atm) where hydrogen-oxygen and hydrogen-air mix- 
tures invariably ignite spontaneously (in other words, at temperatures above 
the junction of the second and third explosion limits). Thus, he tacitly as- 
sumes that there is a "memory" of the second limit - a demarcation between re- 
gions of long and short delay. 

To examine this notion further, ignition delays have been calculated for 
conditions around the lean limit of detonability for hydrogen-air mixtures. 
Chapman -Jouguet Mach numbers were computed as described by Zeleznik and Gordon 
(ref. 18), and then conditions immediately behind the shock were deduced by us- 
ing ideal shock relations for a gas of r = 1.4 (ref. 19). For initial condi- 
tions of room temperature and atmospheric pressure with hydrogen concentrations 
of 12 to 16 percent in air, the Von Neumann spike pressures and temperatures 
are in the ranges 16 to 20 atmospheres and 1100° to 1300° K, respectively. De- 
lays have been computed from equation (30) with rate constants for the initia- 
tion reaction and reactions I, II, and III from reference 4. Rate constants 
for reaction VI were obtained from equation (26) and the constants for reaction 
XI were taken as 


kn = 5. 4x10® 


exp 


-24,000\ 

V RT I 


liters/(mole) (sec) 


which is simply 4.5 times the value suggested by Voevodsky (ref. 12). 
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Results are shown in figure 4. Clearly there is an abrupt decrease in 
ignition lag arouiid a hydrogen mole fraction of 0.14; the ignition delay de- 
creases by two orders of magnitude as the hydrogen concentration is increased 
0.01. This provides further support to Belles’ postulate that the condition 
2kg > kgc^ must be satisfied for stable detonations in the hydrogen- oxygen 

Experimentally, the lean 
limit of detonability of hydrogen- 
air mixtures lies in the range 
14 to 15 percent hydrogen by vol- 
ume (ref. 20) . 


CONCLUDING REMARKS 

It is hoped that the examples 
discussed in the preceding section 
illustrate the usefulness of these 
analytic solutions as applied to a 
variety of problems. Further 
analysis of this sort may help -un- 
cover more detailed information on 
the hydrogen-oxygen reaction. 

Some other applications suggest 
themselves ; 

(l) The concentrations of 
labile intermediates - H, OH, 0, 
H0 2 , etc. - can be computed as a 
function of time throughout the 
induction period. Thus, one can explore the feasibility of detecting these 
species experimentally, for instance, spectoscopically, or with electron para- 
magnetic resonance techniques, or the mass spectrometer. 

(2) By analyzing experimental ignition delay data obtained under a variety 
of conditions, it should be possible to glean information as to the mechanism 
of initiation, that is, the process by which atoms and free radicals are first 
generated from the reactants (by analysis of the term B in eq. (24)). 

(3) By computing delays near the low-pressure limits of detonability, it 
may be possible to gain insight as to the degree of coupling between shock 
front and reaction zone required in a stable detonation. 

(4) If need be, other chemical reactions may be added 
(This may be desirable particularly in the high-pressure - 
region.) There is a restriction, however. Reactions must 
the intermediates or products. Thus, reactions such as H 
cannot be considered. 


to the scheme, 
low- temperature 
be first order in 
+ H 2 0 2 ^H 2 0 + OH 


system. 
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Figure 4. - Ignition delays for Von Neumann spike condition 
in Chapman-Jouguet detonation of hydrogen-air in neigh- 
borhood of lean limit of detonability. 
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(5) Finally, the general method of attack used here should he useful in 
obtaining solutions to the kinetics for other isothermal branched chain explo- 
sions. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, March 26, 1964 
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